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Abstract 

Front propagation in a random environment is studied close to the depinning thresh- 
old. At zero temperature we show that the depinning force distribution exhibits a 
universal behavior. This property is used to estimate the velocity of the front at very 
low temperature. We obtain a Arrhenius behavior with a prefactor depending on the 
temperature as a power law. These results are supported by numerical simulations. 

1 Introduction 

Front propagation in a random environment has in recent years become a generic problem 
describing the motion of interfaces as different as magnetic domain wallsjffl], wetting|| or 
fracture fronts ||. In all cases the balance between a roughening effect due to the quenched 
random environment and the smoothing effect of the elastic interactions leads to a rich phe- 
nomenology. At zero temperature one observes a depinning transition: a threshold force is 
required for unlimited motion. When the front is driven at a force close to the threshold, 
one observes spatio-temporal organization, the front roughness exhibits self-affine correla- 
tions and the critical exponents characterising the propagation of the front can be shown 
to depend on the nature of the elastic interactions [|J. At finite temperature the dynamics 
is controlled by thermal activation and a creep regime characterized by a stretched expo- 
nential dependence of the velocity on the driving force has been predicted by scaling|| or 
renormalisation || arguments. This creep regime has been experimentally observed in weak 
pinning conditions for the propagation of a magnetic domain wall |l| . 

In the following we consider the case of strong pinning, we first give results obtained in 
the framework of an extremal dynamics focusing on the depinning force distribution in the 
vicinity of the macroscopic depinning threshold. We then consider the thermally activated 
propagation of a front driven at a constant force and we show that at low temperature the 
knowledge of the depinning force distribution leads to a simple Arrhenius regime with a 
prefactor that depends on the temperature as a power law. To illustrate these results we 



present in the following numerical simulations performed with two different elastic inter- 
actions: Laplacian (corresponding to the motion of magnetics walls or fluid invasion in a 
Hele-Shaw cell) and long range elasticity (which corresponds to the propagation of a wetting 
or a fracture front). 



2 Force distribution in the extremal model 

We consider the propagation of an elastic ID front through a 2D random environment. The 
position of the front at abscissa Xj is given by hi = h(xj). The random environment consists 
in traps of depth 7$ which are randomly distributed in the direction of propagation. In the 
following we restrict ourselves to the case of an overdamped dynamics in a strong pinning 
situation. When submitted to an external driving force f ex t the depinning criterion for the 
site i can be written as: 



fext + fel[Xi,hi] > 7i (1) 

where the elastic component can be Laplacian or long range. In the particular case of an 
interfacial fracture front, the long range elastic kernel comes from a first order expansion of 
the stress intensity factor || and with periodic boundary conditions it can be written: 

h- — h- 

fei[xi, hi] = . r r > ( 2 ) 

jjH sin [7r{Xi - xrf/L] 

where L is the size of the system. In the Laplacian case, we use a discrete form of the Laplace 
operator: 

fLa P [xi, hi] = - (hi+i - 2hi + hi-i) (3) 

In the extremal dynamics the external force is adapted at every step to the value f ext = 
f c (t) such that only one site can depin. This corresponds to selecting the weakest site io of 
the front: 

fc = lio - fei [x io ,h io ] = mini (7; - f el [x h hi)) . (4) 

The weakest site is then advanced up to the next trap, the elastic forces are updated and 
the process is iterated on the weakest site of the new configuration. The maximum of these 
depinning forces f c over time corresponds to the macroscopic depinning threshold /*. 

The front then presents a multiscale roughness: it remains statistically invariant under 
the anisotropic scaling transformations x — ► Ax, y — ► X^y. The front is said to be self-afnne 
and the exponent ( is called the roughness exponent. This scaling invariance is responsible 
for long range correlations of the height differences on the front. In particular when measured 
over a size d, the height standard deviation grows as a power law: y/< Ay 2 >(d) oc <f>. A 
direct consequence of this self-affinity property is that the height power spectrum is a power 
law of exponent —1 — 2£ where ( ~ 0.35 and ( ~ 1.25 in the fracture front case and in the 
Laplacian case respectively. 

It can be shown that the behavior of the ditribution Q(f c ) in the vicinity of the threshold 
/* exhibits universal features. In order to illustrate this point let us focus on a sequence of 
successive depinning events on the front. Such sequences are known to be time and space 
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correlated and present an avalanche-like behavior Q]. The distance d between two successive 
events is distributed according to p(d) oc d~ a where a = 3 in the Laplacian case and a = 2, the 
exponent of the elastic kernel, in the case of long range elastic interactions. Large separations 
between two succesive events correspond to the jump from one avalanche to another one, the 
larger the jump, the closer the depinning force to the threshold /*. Let us consider the force 
distribution qd(f) conditionned to the size d of the jump. From the knowledge of the elastic 
interactions and of the front roughness statistics we can give estimates of the moments of 
these distributions. The typical height fluctuations on a distance d scaling as <f% the force 
fluctuations scale as d~ b where 6=1 — C — 0-65 in the case of long range interactions or 
b = 2 — ( ~ 0.75 in the Laplacian case. As d increases, qd{f) becomes narrower and closer 
to the threshold /*. In particular we get the linear relationship: 

(f c ) d = f*- Aa d (f c ) , (5) 

where {f c )d and <Jd(fc) are the mean and the standard deviation of the depinning forces 
conditionned to a jump of size d. This relationship allows a precise determination of the 
threshold /*. Moreover all distributions qd are identical up to a rescaling transformation: 

q d (/* - fc) = dfy [(/* - fc)d b ] . (6) 

The knowledge of the distribution p(d) oc d~ a of the distances between successive active sites 
allows to express the depinning force distribution close to the threshold: 

Q(f* -fc)= f x b ~ a 4> \(f* ~ fc)x b ] dx oc (/* - hY , " = a ~ h ~ l , (7) 



vith the numerical values v ~ 0.5 in the fracture front case and v ~ 1.7 in the Laplacian 
case. 

In Fig. 1 we present numerical results obtained in the case of long range elastic interac- 
tions, we see that the low tail contribution corresponds to activity jumps of size or 1. For 
contributions corresponding to larger jumps, the details of the random trap distribution are 
washed out and we can check that the force distributions obey the scaling proposed in Eq. 

Note in addition that the universal behavior proposed in Eq. (|7D only applies in a very 
close vicinity of the threshold /*. The size of the sytem is L = 1024, the trap depth and 
the distance between two traps in the propagation direction are uniformly distributed in the 
range [0 — 1]. The results are averaged over 5 millions iteration steps. A linear fit using Eq. 
(§) gives /* = 0.373 ± 0.001 and we get for the roughness exponent ( = 0.34 ± 0.01. 



3 From extremal to activated dynamics 

Instead of adapting the driving force at every step so that only one site can depin at a time 
we consider now the front propagation under a constant external driving force f ex f At zero 
temperature the front moves freely if f ext > /*, while it will only advance a finite distance 
before being pinned if f ext < /*. This travelling distance diverges when f ext approaches the 
threshold /*. At finite temperature we now allow a pinned site i to depin at any time with 
an Arrhenius probability pi. 
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Figure 1: Left: Distribution of depinning forces (bold) and contributions conditionned by the 
distances between consecutive active sites. Apart from the first distribution (continuous line) 
associated with very short jumps, which is sensitive to the details of the random threshold 
distribution, the other ones present the same tendency: the larger the jump, the closer the 
mean force to the threshold and the narrower the distribution. Right: these distributions are 
scaling invariant, after rescaling according to Eq. @ all these distribution fall on the same 
master curve. These results have been obtained in the case of long range elastic interactions. 



p l = , Aj = f ext + ff- 7i< 0, (8) 

and (assuming all other sites to be frozen) the probability that the site has remained pinned 
during the time r and advances in the interval [r, r + dr] can be written 

1 -x. 1 
PArjdr — —e T ^dr , n = — = e t . (9) 
r % Pi 

If we consider now all sites of the front, the typical waiting time for the first depinning event 
on the front is r* such that: 

^ = E7 = E«*. do) 

' i T * i 

the probability that the depinning takes place at site i being Pi/ Y,Pi = r * / r i- At very low 
temperature the waiting time r* is dominated by the waiting time of the weakest site io 



± = 2_l + V e ^|. (n; 



rn a x 



In the following we focus on this situation of very low temperatures. We consider the 
case of an elastic front driven at a constant force f ext . The (very low) temperature allows 
to introduce a physical time in the extremal model which is adapted as follows. When 
the driving force f ext exceeds the local threshold f c , all sites which fulfill the depinning 
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Figure 2: Velocity of the front after rescaling by the Arrhenius term exp[(/* — f e xt)/T] 
for different driving forces in the case of long range elastic interactions (left) in Laplacian 
elasticity (right). The symbols correspond to the numerical simulations and the lines to 
the expected results from Eq. ( |T2"| with an exponent v = 0.25 (long range) and v = 1.66 
(Laplacian) 

criterion are allowed to advance at the same time. Then the elastic forces are updated to 
take into account the change of front conformation. If however the site reaches a blocked 
configuration where no site is able to depin at zero temperature, the weakest site is allowed to 
advance by thermal activation. The waiting time to is chosen according to the distribution 
(1/Tj ) exp(— to/Ti ) corresponding to the extremal site and the elastic forces are updated 
before the next iteration. 

The distribution of the depinning forces being known in the vicinity of the macroscopic 
threshold /*, we can estimate the creep velocity at very low temperature via the average 
waiting time of an activated step: 
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so that we expect at low temperatures: 



,{f exU T)=T-^ 



cxp 



/ f ext 



T 



(13) 



and a simple Arrhenius behavior when T 3> /* — f ext . 

In Fig. 2 we present the numerical results obtained for a front driven at different forces 
below the threshold. The symbols correspond to the velocity results after correction by the 
Arrhenius term and the lines to the expected behavior following Eq. [12]. For the latter 
expressions we used the values v = 0.25 in the long range case and v = 1.66 and Laplacian 
case to compare with the expected values v = 0.5 and v = 1.7 respectively. These differences 
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may be due to the fact that the universal behavior of the depinning force distribution is only 
valid in a region close to the threshold (Af < 0.1 in our case, see Fig. 1) and for lower 
driving forces the system becomes more sensitive to the details of the distribution. 

This behavior is different from the stretched exponential described in |5], |6], [7], ||; the main 
reason is probably that although we consider very low temperatures, our study is restricted 
to the very close vicinity of the threshold while the usual scaling obtained for creep motion 
corresponds to driving forces far below the threshold. Note in addition that at very low 
temperature, the time scales of depinning become very large and that it may be necessary 
to take other mechanisms into account. This is for example the case for crack propagation 
in glass for which the subcritical motion of the front competes with a stress assisted ion 



interdiffusion phenomenon which tends to locally reinforce the material ||10|| . Within these 
restrictions, our knowledge of the depinning force distribution gives us the evolution of the 
front velocity with temperature by way of the front roughness, thus relating a dynamical 
behavior to a geometrical property. 
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